arXiv:1507.07552v2 [nucl-th] 22 Oct 2015 


INT-PUB-15-034 


Matching pion—nucleon Roy—Steiner equations to chiral perturbation theory 

Martin Hoferichter, 1,2,3 Jacobo Ruiz de Elvira, 4 Bastian Kubis, 4 and Ulf-G. Meifiner 4,5 

1 Institut fur Kernphysik, Technische Universitat Darmstadt, D-6J289 Darmstadt, Germany 
e ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum fur Schwerionenforschung GmbH, D-6f291 Darmstadt, Germany 
3 Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550, USA 
* Helmholtz-Institut fur Strahlen- und Kernphysik (Theorie) and 
Bethe Center for Theoretical Physics, Universitat Bonn, D-53115 Bonn, Germany 
5 Institut fur Kernphysik, Institute for Advanced Simulation, 

Julich Center for Hadron Physics, JARA-HPC, and JARA-FAME, 

Forschungszentrum Julich, D-52425 Julich, Germany 

We match the results for the subthreshold parameters of pion-nucleon scattering obtained from a 
solution of Roy-Steiner equations to chiral perturbation theory up to next-to-next-to-next-to-leading 
order, to extract the pertinent low-energy constants including a comprehensive analysis of systematic 
uncertainties and correlations. We study the convergence of the chiral series by investigating the 
chiral expansion of threshold parameters up to the same order and discuss the role of the A(1232) 
resonance in this context. Results for the low-energy constants are also presented in the counting 
scheme usually applied in chiral nuclear effective field theory, where they serve as crucial input to 
determine the long-range part of the nucleon-nucleon potential as well as three-nucleon forces. 

PACS numbers: 13.75.Gx, 11.55.Fv, 12.39.Fe, 11.30.Rd 

Keywords: Pion-baryon interactions, Dispersion relations, Chiral Lagrangians, Chiral symmetries 


INTRODUCTION 

Chiral symmetry of QCD, the invariance of the QCD 
Lagrangian under chiral rotations of the quark fields in 
the chiral limit of vanishing quark masses, is a power¬ 
ful tool to elucidate the properties of strong interactions 
at low energies, where QCD becomes non-perturbative. 
This chiral symmetry is known to be broken sponta¬ 
neously and explicitly, with the appearance of almost 
massless pseudo-Goldstone bosons, the pions. By ex¬ 
panding systematically around the chiral limit of van¬ 
ishing quark/pion masses, one obtains an expansion in 
momenta and quark masses, with non-analytic terms pre¬ 
dicted and the effects of high-energy physics incorporated 
in low-energy constants (LECs). These LECs appear in 
different physical processes, so that once fixed in one pro¬ 
cess, they can be used to predict others. In particular, 
one can derive low-energy theorems that relate differ¬ 
ent observables, at a given order in the chiral expansion. 
This approach, Chiral Perturbation Theory (ChPT), was 
pioneered in the meson sector in [U-Q , and manifold ex¬ 
tensions have been worked out over the last decades. In 
particular, it has been extended to the single-baryon sec¬ 
tor, see [l-[ll, with pion-nucleon (irN) scattering as one 
of the most fundamental applications SB®. 

However, as first pointed out in [l6H18l |. constraints 
from chiral symmetry are by no means limited to sys¬ 
tems with at most one nucleon: once so-called nucleon- 
nucleon (NN) reducible contributions are separated, the 
remaining irreducible parts of the NN potential again 
permit a chiral expansion, despite the non-perturbative 
nature of the NN interactions. In this way, Nuclear 
Chiral Effective Field Theory (ChEFT), the extension 
of ChPT to the multi-nucleon sector, has been developed 


as a powerful tool for a systematic model-independent 
approach to nuclear forces, see 19j,l20j for recent reviews. 


One particularly valuable feature of ChEFT concerns the 
prediction of a hierarchy between two- and multi-nucleon 
forces, with the NN interactions starting at leading order 
(LO), three-nucleon forces are predicted to enter at next- 
to-next-to-leading (N 2 LO) order, and even higher forces 
are accordingly suppressed (2lj| . Further, the LECs ap¬ 
pearing in the expansion relate different processes. In 
fact, the LECs that appear in nN scattering determine 
the long-range part of the NN potential and the three- 
nucleon force. Accordingly, if sufficiently precise informa¬ 
tion on nN scattering were available, the required input 
could be immediately used in multi-nucleon applications. 
This is the aim of the present Letter. 

Such improved input for the nN LECs is becoming 
increasingly urgent, since with higher orders in ChEFT 
being worked out the uncertainties in the nN LECs are 
starting to significantly contribute to the error budget in 
some observables, see e.g. 22j. In the past, several strate¬ 
gies have been pursued: extractions from nN scattering 
data, either in terms of phase shifts H, Hf or cross sec¬ 


tions M , determinations from NN observables [H, I_ , 

or a combination of both 271. Moreover, in Q the match¬ 
ing with a reconstructed dispersive nN amplitude was 
performed in the subthreshold region where ChPT is ex¬ 
pected to converge best, but the extrapolation from the 
physical region still required input from nN data (sim¬ 
ilarly, while starting from the subthreshold region, the 
LECs are determined from fits to phase shifts in Il3|). 
Given that the long-range contributions are entirely de¬ 
termined by nN physics, nN scattering provides the 
cleanest access and offers, at least for most LECs, also 
the highest sensitivity for their extraction. However, such 
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a program has been hampered by inconsistencies in the 
low-energy 7 tN data base, as exemplified by contradict¬ 
ing partial-wave analyses, the Karlsruhe-Helsinki 28} 29} 
and the GWU/SAID solutions |30|. 

In 7T7T scattering, a similar situation prevailed until the 
consequent use of Roy equations [31], a combination of 
constraints from analyticity, unitarity, and crossing sym¬ 
metry in the form of coupled integral equations for the 
partial waves. This significantly advanced the knowl¬ 
edge of the low-energy nn phase shifts [32, 331. Indeed, 
the matching to ChPT then allowed for a very precise 
determination of the pertinent 7T7 t LECs [34J. Mean¬ 
while, Roy-equation techniques have been extended to 
other processes 35, 36], in particular, a similar program 
has been pursued for irN scattering based on Roy-Steiner 
(RS) equations 37 - 4(| , making use of a high-accuracy ex¬ 
traction of the 7 tN scatterin g lengt hs from pionic atoms 
as an additional constraint [4li445j| . In this Letter, we 
work out the consequences of our RS solution for the 
nN LECs by matching the RS and the ChPT represen¬ 
tation of the nN amplitude in the subthreshold region. 
The main advantages of such an approach are the follow¬ 
ing: first, the nN amplitude in the subthreshold region 
is a polynomial in the Mandelstam variables (apart from 
the Born terms), so that the chiral series is expected to 
converge best there. In contrast to 0, we do not need 
additional input from the physical region, as in our case 
the subthreshold parameters follow from the RS solution 
alone. Second, the matching amounts to equating the 
subthreshold parameters from jitl 06] with their chiral 
expansion, which reduces the determination of the LECs 
to an algebraic problem. Third, we can use the compre¬ 
hensive error analysis performed in 4C], 461, which trans¬ 
lates to a full covariance matrix for the extracted LECs. 


SUBTHRESHOLD PARAMETERS 
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TABLE I: Subthreshold parameters from the RS analysis [4(il , 
0 . 


index are understood to be related to the 7r ± p —> tt± p 
charge channels according to 


X^r 2 {X,. 


p—t7V p X 7r +p_). 7r -rp 


)> 


(4) 


for A' e {D, B ,..and the nucleon and pion mass are 
identified with the masses of the proton and the charged 
pion, respectively, see and {rMiO} for a discussion of 
the pertinent isospin-breaking corrections. As mentioned 
above, once the Born terms are subtracted, the amplitude 
in the subthreshold region becomes a polynomial in v and 
t. A particularly convenient representation is provided by 
the subthreshold expansion 


D ± (v,t) = ( 1 ) f; d±y m t n -, 

' ' n,m =0 

B ± (u,t)= b ^y m t n , (5) 

'■ ' n,m=0 

where the upper/lower entry corresponds to I = ±, and 
the Born-term-subtracted amplitudes are defined as 

X±{v,t) = X^t) - X±(v,t), X e {D,B}, (6) 


We start by specifying conventions for the process with 


ir a (q) + N(p) ^ ir b (q') + N(p), (1) 

with pion isospin labels a, b and Mandelstam variables 

s = {p + q) 2 , t = (p'-p) 2 , u = (p — q') 2 , (2) 

fulfilling s + t + u = 2 m 2 N + 2 M 2 . We parameterize the 
scattering amplitude as 

T b >,f) = 6 ba T+(v,t) + i[r 6 ,r a ]T-(M), 

T J (i/,t) = u(p')|n 7 (ci) - ^^■B I (u,t)^u{p), (3) 

where v = (s — u) / (Attun), the isospin index I = ± refers 
to isoscalar/isovector amplitudes, win and M n to the nu¬ 
cleon and pion mass, and r a denotes isospin Pauli matri¬ 
ces. Throughout, the amplitudes with a definite I = ± 


b p- y,t) = s 2 ( 
D py,t) = — 


1 1 A 

—2--F —2 - 

m N — s m N — u J 

■ ^ +vB± v (y,t), 




(7) 


where g denotes the 7 tN coupling constant. 

For the matching to ChPT at N 3 LO (complete one- 
loop order) we need the 13 subthreshold parameters listed 
in Table H The solution of the RS equations is obtained 
by minimizing a y 2 -like function, defined as the difference 
between left- and right-hand side of the equations on a 
grid of points, with respect to the subtraction constants 
and the low-energy phase shifts. Most of the subthresh¬ 
old parameters listed in Table [T] already appear as sub¬ 
traction constants of the RS system, and thus follow as 
output from the RS solution, while the remaining ones, 
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, . and o?o 2 ! are calculated from sum rules after¬ 


wards. The uncertainty estimates include a number of 
effects: first, the RS equations are valid only in a finite 
energy range below the so-called matching point and only 
a finite number of partial waves are included explicitly in 
the solution. We varied the input for the matching condi¬ 
tion as well as for the energy region above the matching 
point and higher partial waves, both regarding different 
partial-wave analyses and truncations of the partial-wave 
expansion. Furthermore, we varied the input for the irN 
coupling constant within g 2 /( 47 r) = 13.7(2) [44, 451 and 
investigated the sensitivity to the parameterization of the 
low-energy phase shifts used in the solution. Second, 
we observed that the RS equations are more sensitive to 
some subthreshold parameters than others. To account 
for this effect, we generated a set of solutions correspond¬ 
ing to different starting values of the y -minimization, 
while imposing sum rules for the higher subthreshold pa¬ 
rameters, and took the observed distribution as an addi¬ 
tional source of uncertainty. Third, we propagated the 
errors in the scattering lengths, which crucially enter as 
constraints in the minimization, to the results for the 
subthreshold parameters. Taking everything together we 
obtain a 13 x 13 covariance matrix that encodes uncer¬ 
tainties and correlations of the 13 subthreshold parame¬ 
ters relevant for the matching to ClrPT. 


CHIRAL EXPANSION 


The chiral expansion for the subthreshold parameters 
is spelled out explicitly in 0 , in particular 


,+ 2M^(2ci-c 3 ) ffj(3 + 8ffj)Mg 

00 F 2 + 64ttF4 

, f4 f 16eu 3g^(l + 6g\) _ 2ci - c 3 \ 

71 1 F 2 647r 2 F4 TOj v 16t r 2 F;) J ’ 

1 , 4 M 2 (di + d 2 + 2 d 5 ) 

vo - 2F 2 + F 2 + 48 n 2 F 4 

_ M 3 ( 8 + 12g^ + llg^ _ 4ci +g 2 A {c 3 - c 4 ) 
71 | 1287T F^rriN 47t F£ 


( 8 ) 


where Cj, J,;, and e, denote the NLO, N 2 LO, N 3 LO nN 
LECs, respectively, F' n the pion decay constant, and gA 
the axial coupling of the nucleon. The conventions for 
the &i correspond to the general classification 0 and 
the a have been redefined to absorb a quark-mass renor¬ 
malization, see 010 . Finally, the expressions in <[ 8 |) 
follow the standard counting in the single-nucleon sec¬ 
tor, where the expansion parameter is given by 0{jp) = 
{p, Mtt}/A b, for momenta p and the breakdown-scale 
Ab ~ A x ~ 47 tF„ ~ tun ~ M p ~ 1 GeV. In contrast, 
the breakdown-scale in few-nucleon applications is typ¬ 
ically lower, Ab ~ 0.6 GeV, so that relativistic correc¬ 
tions are often counted as {p, M n }/mN = 0(p 2 ) 53- 
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TABLE II: Results for the nN LECs at NLO, N 2 LO, and 
N 3 LO (standard and NN counting only differ at N 3 LO, 
except for NLO in c 4 , which in the NN scheme becomes 
2.44(3)). The results for the a , di, and e; are given in units 
of GeV -1 , GeV -2 , and GeV -3 , respectively. 


As a consequence, in this counting one would drop the 
1/mjv suppressed terms in | 8 ]). In this Letter, we con¬ 
sider both counting schemes, which we will refer to as 
standard and NN counting in the following. The full set 
of subthreshold parameters can be easily inverted for the 
LECs, for the result at different chiral orders see Table ITT1 
(masses, F n . and gA are taken from 0). The errors as 
propagated from the subthreshold parameters are tiny 
compared to the shifts observed between chiral orders: 
clearly, the dominant uncertainty now resides in the chi¬ 
ral expansion. For completeness, we also quote the N 3 LO 
correlation coefficients, see Table HTT1 Note that this ta¬ 
ble contains the correlation matrices for the standard and 
the NN counting and therefore appears asymmetric. 


In general, the values for the LECs are expected to be 
0 ( 1 ), e.g. Ci ~ qa /Ab 0 and significant departures in¬ 
dicate the presence of additional degrees of freedom. In 
the case of C 2-4 the main origin of their enhancement is 
well understood: while other resonances do contribute as 
well, it is primarily the presence of the A(1232) resonance 
that makes these LECs take unnaturally large values 52 - 


54J. Following [54], we extract the A contributions to the 


individual subthreshold parameters from the correspond¬ 
ing tree-level A-exchange diagrams and convert the re¬ 
sult to the LECs. For the numerical analysis we use the 
A coupling constant g^NA = 1-2 [0, which lies right 
in the middle of the range 1.05 (extracted from the A 
width [55J) and 3<7 a/( 2\/2) = 1-35 (predicted by large 
N c 0). Keeping the full pion-mass dependence, we ob¬ 
tain the values shown in the first column of Table IIVI 
while the second column corresponds to the leading ex¬ 
pansion in M,r and ijia — niN- Only in the latter case 
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TABLE III: Correlation coefficients at N 3 LO in standard (upper-right triangle) and NN (lower-left triangle) counting. 
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TABLE IV: A contributions to the nN LECs, in GeV units, 
for the full A-exchange diagrams (first column) and to leading 
order in M n and mA — mw (second column). 


one recovers the relation = — cf = 2cf 153 . 

While the A can indeed explain a significant por¬ 
tion of the physical values of the Cj, its effect is too 
small to explain the large numbers for the di that ap¬ 
pear at N 3 LO (except for ds). The origin of this large 
shift can be traced back to the terms proportional to 
9a(°3 ~ c 4 ) ~ —16 GeV -1 in dg 0 in © (and similarly 
in d/g, dy 1: and 6 q" 0 ). These terms mimic loop diagrams 
with A degrees of freedom. Our results show that if the 
A is not included explicitly, such contributions lead to 
a substantial renormalization of the LECs. Indeed, if 
we drop the C 3 — C 4 loop terms, the di are reduced to 
^ 1 + 2 , 3 , 14-15 = (2.2,-3.9,-2.6) GeV -2 , in good agree¬ 
ment with the expectations from Table IIVI 


THRESHOLD PARAMETERS 


With the LECs determined by matching to the sub¬ 
threshold expansion, it is important to check how well 
the chiral series converges in other kinematic regions. 
A prime test case is provided by the S'-wave scattering 
lengths aij ), : they are known very precisely from pionic 
atoms 44|, [45| ■ In the isospin conventions © their values 


are a[} + = —0.9(1.4) and ag + = 85.4(9) (always in units 
of 10 -3 M -1 ), and the problem is still purely algebraic. 
The fourth-order expressions for their chiral expansion 
were first given in [5 71 ] , in our conventions they read 

+ _ M l f 3glm N M„ g\Ml 
a ° + ATtF%(m N + M^X 647 tF% 16 m 2 N 

- j [9a + 8toat(2ci - c 2 - c 3 )] 

+ M 2 [ - 16cic 2 + di 8 g A + 16m N (eu + ei 5 + ei 6 )] 

_ Ml [8 - 3 g\ + 2 g\ + 4mjy (2ci - c 3 )] 

64t t 2 F% 

m^M n f Ml g\Ml 

a °+ = 8nFZ(m N + M^l 1 + 8 + 4 ^“ 

+ 8M 2 (di + di + d 3 + 2ds)^. (9) 

Fixing the only new LEC, dis, from the Goldberger- 
Treiman discrepancy, 

Jls = f 1 - J ^ L ~] = -°- 44 ( 24 ) GeV-2 > ( 10 ) 

2M- \ rnwgA ) 
we obtain the following results 

a+ + = {-23.8, 0.2, -7.9}, {-14.2,0.2, -1.4}, 
a~ + = {79.4,92.9,59.4}, {79.4,92.2,69.2}, (11) 

where the first/second array refers to the standard /NN 
counting and the three entries to NLO, N 2 LO, N 3 LO. It 
is not surprising that the chiral expansion in the isoscalar 
sector is slow, after all its LO vanishes. Unexpectedly, a 
similarly slow convergence is also found for a^ + , whose 
low-energy theorem at LO is tantalizingly close to the full 
answer, while the agreement in both counting schemes 
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deteriorates when going to fourth order. The largest 
part of this discrepancy can be attributed to the A loops 
discussed above, i.e. for a^ + the largest portion of the 
C 3 — C 4 terms does cancel between d\ + d 2 and ds, but the 
cancellation is incomplete and the remainder spoils the 
agreement with the pionic-atom value. 

This example shows that in a theory without explicit 
A degrees of freedom the LECs determined in a partic¬ 
ular kinematic region do not necessarily ensure conver¬ 
gence in the full low-energy domain. However, especially 
when going to higher orders, including the A explicitly 
becomes extremely challenging, so that in practice the 
A-less approach can be pushed to higher orders and it re¬ 
mains to be seen if in the end the A-full or A-less theory 
proves more efficient. We argue here that for A-less ap¬ 
plications in NN scattering matching at the subthreshold 
point is the preferred choice: the two-pion-exchange dia¬ 
grams can be reconstructed from nN scattering by means 
of Cutkosky rules 0] (see 0 for recent applica¬ 
tions of this approach), with spectral functions involving 
irN amplitudes either directly evaluated at or weighted 
towards zero pion center-of-mass momenta [ 0 , which 
translates to s = m 1 2 3 N — M 2 . Moreover, for physical val¬ 
ues of the momentum transfer t in NN scattering the 
Cauchy kernels in the spectral integrals become largest 
for t = 0. Since the corresponding combination of (s, t) is 
much closer to subthreshold (m%+M 2 , 0 ) than threshold 
((mjv + M n ) 2 , 0) kinematics, we conclude that the LECs 
to be applied in nuclear forces should be extracted from 
the subthreshold point in 7 tN scattering instead of the 
physical region. In the present Letter we have presented 
such an extraction based on a comprehensive analysis of 
low-energy ttN scattering in the framework of RS equa¬ 
tions. The corresponding LECs clearly defined at a given 
chiral order will be valuable for assessing the uncertain¬ 
ties from the long-ran ge p art of the nuclear force in future 
ChEFT calculations 
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